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, ^ . Abstract 
\D 

\ The renormalization group equations(RGEs) of non-universal soft supersymmetric 

QQ ' breaking terms with CP violating phases are analyzed in this paper. We obtain the 

. analytic solutions of RGEs by directly solving the RGEs themselves. Compared with the 

method of spurion expansion our approach proves to be simple and succinct, and easy to 
extend to the case of complex parameters. With the analytical forms of the solutions we 
obtained the infrared quasi fixed point behavior of soft terms are analyzed and it turns 
out to support the notion in scenarios with CP violating phases. 
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1 Introduction 



As an excellent candidate to embed the standard model(SM) in a more fundamental theory, 
the minimal supersymmetric standard model(MSSM) has many good features. It ingeniously 
tackles the abominable gauge hierarchy problem from which ordinary unification theories suffer 
[|^, provides a mechanism that supergravity models all share which breaks the electroweak(EW) 
symmetry dynamically via radiative corrections, naturally provides the lightest supersymmetric 
particle(LSP) as a candidate of the dark matter in astrophysics and cosmology^, etc. In 
spite of these theoretical virtues, however, the MSSM suffers from great uncertainty which 
arises from the large number of free parameters describing the soft supersymmetry (SUSY) 
breaking. One may greatly reduce the parameter space and make the theory much more 
predictive via adopting the universality conditions at high energy scale, e.g., one may choose the 
minimal supergravity model (MSUGRA). Furthermore, with the aid of the concept of infrared 
quasi fixed point (IRQFP), further reduction of parameter space has been observed in many 
works ^ 1^, 0. It reveals the screening effect for soft breaking parameters implied entirely 
by the renormalization group equations(RGEs) and the large top Yukawa coupling, and with 
this one may rely on two0 or even one^ free soft SUSY breaking terms, i.e., Mq and M1/2 
or only Mi/2- In these cases one can work in models of highly predictive power leaving some 
of the universality boundary assumptions. This IRQFP analysis was extended to the large 
tan/3 scenario]^ in which the whole set of Yukawa coupling of the third generation has to be 
considered. The result shows that the theory also exhibit the IRQFP for large tan/? and allows 
one to predict SUSY mass spectra as functions only of M1/2, but the analysis turns out to be 
much more complicated. 

While hindered in this situation, a crucial observation was made in Ref. p] which showed 
the intrinsic connection between the solutions for the RGEs of soft terms and those of the gauge 
coupling constants and Yukawa couplings. It tells that one can get the solutions of RGEs for 
the soft SUSY breaking terms via substituting the modified expressions for gauge and Yukawa 
couplings in the solutions of the RGEs for gauge and Yukawa couplings P, P|, pil[]. 

«. =^ ai = a,{l + Mir] + Mif] + 2MiMirif]), t] = 9'^, f] = 9'^, 
n n = n(l - AkV - Akf] + AkAkVV + ^kVV) , 

and expanding over the Grassmannian parameters 6 and 6, which is exactly the spurion field 
description of soft SUSY breaking The technique was immediately used to probe the 



implication of RGEs for the soft SUSY breaking terms in large tan/3 scenario, but one faced 
with the problem of finding good analytical forms for Yukawa couplings. This problem was 
solved in Ref. [|12| in which an integral form for Yukawa coupling evolution was given which was 



shown to be a convergent scheme for numerical calculations and a convenient and concise form 
for qualitative analysis. In view of this solution, the general analytical forms for the solutions 
of the soft SUSY breaking terms were obtained in this complicated case via the spurion field 
expansion and the IRQFP feature was analyzed for non-universal boundaries of soft terms |]T3| . 



Although the solution for soft SUSY breaking terms look concise and simple, the procedure 
described above, especially the actual calculation of spurion expansion is still complicated due 
to the a little bit complicated iterative solutions for Yukawa couplings |12| . So in Ref. [jl3| the 



author restrict their attention to real trilinear couplings and real gaugino masses. However, 
in a lot of works on superstring theory and M-theory phenomenology non-universal complex 
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soft terms appear naturally which provide new sources of CP violation. And the SUSY 
CP violating phases may not be constrained to be small due to the cancellation mechanism of 
SUSY contributions to electron and neutron electric dipole moments (EDME and EDMN), and 
one has extensively discussed the phenomenological implications of large phases of /i, gaugino 
masses and trilinear terms while satisfying the EDME and EDMN constraints [0. Therefore, 
it is important to extend the study of non-universality of soft terms to the complex parameter 
case. In this paper we will show that the convenient iterative forms[jl2|, |13[|of the solutions of 
RGEs can be obtained very succinctly in the approach of directly solving RGEs and extend to 
the complex non-universal soft terms. Our results show that the real and imaginary parts of 
trilinear couplings have infrared quasi fixed points respectively if the initial values of Yukawa 
couplings are large enough ( at least > a^). 

The paper is organized as follows. In section two we will show briefly how one can directly 
solve the RGEs and arrive at the forms in Ref. ||l3l, and extend the analysis to more general 
case, complex parameter case. In section three we will analyze the IRQFP in the case with 
SUSY CP violating phases and discuss the phenomenological implications. Finally, conclusions 
are drawn in section four. 



2 Solution of the RGEs 

The one loop RGEs for gauge and Yukawa coupling and soft breaking terms can be written 



as 

di = -bia^, (1) 
Yk = Yk{J2'^kiai-^"'kiYi), (2) 

i I 

Mi = -haiMi, (3) 
= -Y.(^kiaiMi-Y,akiYiAu (4) 

i I 

tk = 2^Cfc,«,|M,|2-^afc,l^(S, + |A;|2), (5) 

i I 

where ctj = gf/lGir'^, i=l,2,3, Yi = yf/lGir'^, l=t,b,r, ■ = d/dt, t = log Mqjjj^/Q'^ and 

v-i ~ 2 I ~ 2 I 2 

= rhl^ + mis + mjj^, 

= mis + mis + mjj^, 

hi = {33/5,1,-3}, 

Cti = {13/15,3,16/3}, Cb, = {7/15,3,16/3}, c^^ = {9/5, 3, 0}, 

au = {6,1,0}, aM = {1,6,1}, a^; = {0,3,4}. 

The solutions for a^'s and Mj's can be easily obtained to be: 

«i = Mi = (6) 



2 



For Yfc's, solutions turn out to befT^, [13 



1 + akkYk lo Uk 



Yy = (7) 



where UkS are defined iteratively as 

Uk = Ekllil + auYi^ J^i)-'""^''", (8) 



with 



Ek = f[{l + ha'^tr'='/'\ (9) 

i=l 

To find the analytic solutions of the renormalization group equations for AkS and S^'s, we 
rewrite them as the form of the Bernoulli equation. Consider the trilinear coupling Ak first. 
We can rewrite Eq. (4) to be 

^ = - akkYkAk, (10) 

with Cfc defined by 



dck 
dt 



Eq. (10) can then be solved in the standard way. The solution is 



CkiaiMi + akiYiAi, 

i ly^k 

ek{t)\t=o = 0. (11) 



A^l + ttkkYk lo Uk^ k 
1 + akkYk /o 



Ak = -ek + ^ — (l^j 



where A^ = Ak\t=o- In order to solve Eq. (11), we rewrite Eq. (10) as 

d{Ak + ek) 



dt 

which means that 



akkYkAk (13) 



VA— ( A -L ^ aO\ _ yO ^k lo Uk lo UkCk , . 

YkAk- {Ak + Ck- Aj.) - Y^—— — ^—^ . (14j 

Ofcfc l + akkYkhUk 

Integrating Eq. (11) over t and inserting Eq. (14) in it, we find the iterative integral equations 
for Cfc's 

ek = t}^ CkiaiMi + 2^ aki ji — . (15) 

i l^k ^1^1 +0'llJoUl 
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For Efc's, the situation is similar but a little more complicated. Prom Eq. (13) we can write 
d\Ak + ek\^ 



dt 



= -akkYkAk{Ak + Cfc)* - akkYkAKAk + Ck), 



or 



d{\Akf + Alek + Akel) _ d\ek\' 



dt 



dt 



akkYklAk]"^ - akkYk{\Akf + Alck + Akel). 



Define 



^fc — Sfc ~ l^fcT ~ -^k^k ~ -^k^kJ 



then we find from Eqs. (5) and (16) that 

dEfe dC,k 



dt 



dt 



— o^kkYk^k 



(16) 



(17) 



(18) 



which is the same as Eq. (10) with Ak and Ck substituted by and ^k respectively. In Eq. 
(18) ^kS are defined by 



d^k 
dt 



2J2ck^ai\Mi\' -J2akiYi{^i + \Ai\^) + 



li^k 



d\ek? 
dt ' 



Then from Eqs. (17) and (18) we get 

Sfc = Cfe + I + ^feCfe + Akel - 



\Al\^-Yl^akkYniukik _ 
1 + akkY^l /o Uk 



(19) 



(20) 



The iterative equations for ^jt's are derived in a way similar to that for Cfc by noting that Eq. 
(18) can be rewritten as 



dEfc _ dik (i\ek\' 



dt dt 



dt 



akkYk^k + 



which leads to 



rn(Efe + \Ak\^) = - — (Efe - E° - Cfe + \ek\^). 

Jo dhh 



'0 Qfcfc 
From the above equation and Eq. (19), we find 

ik = I E ^ki^iM^ r + 2t E I r - E I 1 2 

i i i 

+t 2^ Cfei Qii 2^ aki Ti + t ^ CkiaiM^ ^ — 

i f^k ^/V + aiiloUi T J^k ^/V + auJoUi 

^° /o - So ^1(^1 \ V- ^? /o - /o w/e; 2 



+ 



+ E '^fc^'^'i 



E^w 

l^k 



(I + E°) - A° tx.e; - uid + /o* 

l/r,o + a«/o*M^ 



(21) 
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The Eqs. (12), (15), (20) and (21) are our main results. When Mj and Ai are real they reduce 
to the results given in Ref. [0, as it should be. The Eqs. (12) and (15) are the same as those 
in the case of real Ai^s and Mj's so that it follows that Ai is independent of its initial value when 
the initial values of all three Yukawa couplings are large enough {Yi^ larger than ). That is, 
the real and imagine parts of Ai divided by M3 have the IFQFPs, respectively. By inspecting 
Eqs. (20) and (21), one is led to that possesses an IRQFP if the initial values of all three 
Yukawa couplings are large enough. In conclusion Yukawa couplings Yi and soft terms Ai and 

have IRQFP behaviors in the complex non-universal Ai and Mi case when the initial values 
of Yukawa couplings, yj°'s, are large enough. 

It is worth to note that having the analytic solutions of RGEs, it is easy to find the RG 
invariants. For example, if we define 

i=l 

with E'^ = E't^{t = 0) and < = u^(t = 0), then from Eqs. (7), (12) and (20) we obtain the RG 
invariants 



Elk — 77 Clkk / Mfe, 

Yk Jo 

E2k = (Ak + ek-akk^ [ Uk<^k^ fl - cikk^ [ u'k] , 

\ U)^ Jo /V Mfc Jo J 

Ezk = (^k-Ck-\Ak\'^-Alek-Akel + Akk^[ MfcC^Vl - a^A;^ / m'^) , 
V Uf, Jo / \ Uf. Jo / 

where k = t,b, r, which extend the results in Ref. to the case with multi Yukawa couplings. 

3 Numerical Study of Dependences of Soft Terms on 
Their Initial Values 

In order to see the dependences of soft terms on their initial values one can write 



Akit) = Y.^imUT.d1ml (22) 

I i 

S.(t) = E fam<? + E^?^(t)A°*A° + E hmM^A'}* + c.c.) + E 4(t)M°Mf (23) 

Oi IV i,l i,j 

where the superscript of A^, Mf etc. means that they are the initial values, i.e., the values 
at unification scale, the asterisk, *, means the complex conjugate, g^i = gfq and 1^^ = In 
Eq. (^) a runs over all the third generation scalar quarks and Higgs scalars. The coefficients 
in the equation depend on the initial values of Yukawa and gauge couplings and , of course, are 
determined by the solutions of RGEs, (12), (15), (20) and (21). 
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The IFQFPs exist for any values of tan/3 if the initial values of all three Yukawa couplings, 
F^P's, are large enough (say, larger than the gauge coupling at GUT scale), as shown in the last 
section. However, that whether the condition can be satisfied is in fact dependent of tan/3 if 
we want to make the theory realistic, i.e., to impose the requirement that the third generation 
quark masses are given by experiments. Therefore, for numerical calculations, we choose the 
procedures as follows. First, run gauge couplings from scale up to find the scale where 
gauge couplings unify with the SUSY threshold effects taken at Msusy = 400 Gev. The values 
of gauge couplings at mz are taken directly from Ref . fl^ . For three Yukawa couplings, we take 
corresponding current masses of the third generation quarks from Ref. ||T^ and calculate their 
pole masses, then run them to rriz scale to find the Yukawa couplings at that scale with tan/? 
fixed (hence the SUSY corrections to the pole mass of top are not included in our calculations). 
The values of Yukawa couplings at GUT scale, Y^^s, are found by running their energy scale 
value up by using Eqs. (7,8,9) and again the SUSY threshold are taken into account at Msusy- 
The second step is to run the gauge couplings, Yukawa couplings, gaugino masses, trilinear soft 
SUSY breaking terms and T,i{l = t,b and r) down to the low scale (e.g., the Msusy scale ) 
with the aid of the integral equations (6-9), (12), (15), (20), (21). The masses of the sfermions 
of the third and the Higgs bosons can be got via S;(/ = t,b and r) as shown in Ref. [0. With 
such a procedure the dependences on Yukawa couplings of the coefficients in Eqs. (22p3D are 
translated into the dependence on tan/?. We show below an example of numerical results at 
ts = log{M^uT/MlusY) = 63.28 for tan/3 = 58: 
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0^ 



At{ts) = 0.19118A° - 0.04580v4^ + 0.01059A° - 0.02728M° - 0.20187M° - 1.45580M^ 
Ab{ts) = -0.04019A° + 0.06052A|^ - 0.03687A° - 0.00168M° - 0.14498M° - 1.324111/°, 
A^{ts) = 0.03346A° - 0.20839A° + 0.28668A° - 0.08353M° - 0.22527M° + 0.470137^°, 
St(t5) = -0.04580(m^3)2 - 0.03520(m$^i)2 + 0.19118(m2^2)^ + 0.01059(m2^3)^ + 0.01059(m|3) 
+0.14538(m^3)2 + 0.19118(m^3)2 - 0.15170|A°|2 + O.OOMOl/lj^l^ + 0.00227|A°p 
+0.02344(A°A||* + A°*A||) - 0.00314(4° + A°*A°) - 0.00071 (A^A°* + A|J*A°) 
+0.00841 (M°A°* + M°M°) - 0mU9{M^Al* + M^*A°^) + 0.00000(M°A°* + M^M") 
+0.04673(M°A°* + M°*A°) - 0.00865(M°A°* + M°M°) + 0.00072(M°A°* + M°M°) 
+0.18528(M3°A°* + M^*A^t) - 0.03551 (M^Al^* + M^*Al) + 0.00515(M3°A°* + M^M") 
-0.00418(M{'M°* + M°*M°) - 0.01949(M{'M3°* + M^M") - 0.13369(M°M3°* + M^*M*^ 
+0.03411|M°p + 0.33227|M°P + 4.96704|M3Y, 
Sb(ts) = 0M052{m%f + 0.02364(m^i)2 - 0.04019(m^2)^ - 0.03687(m2^3)^ - 0.03687(m^3)2 
+0.02032(m° 3)2 _ o.04019(m^3)2 + 0.00420|A°|2 - 0.03908|A°|2 - 0.00028|/l^p 

+0.00441 (A°A°* + A°*A°) - 0.00225(A°A°* + A°*A°) + 0.00859(A°A°* + A°*A°) 
+0.00017(M°A°* + M°*A°) - 0.00028(M°A°* + MfA^) + 0.00089(M°A°* + MfA^) 
-0.00243(M2M°* + M°M°) + 0.00974(M°A|^* + M°M|J) - 0.00246(M°A°* + M°*A°) 
-0.01353(M3°A°* + M3°M°) + 0.05688(M3°A|^* + M^M^) - 0.03167(M3°A°* + M^M") 
-0.00146(M{'M°* + M°*M°) - 0. 00279 (MfM^* + MfM^) - 0.09627(M°M3°* + M^*Mr^ 
-0.00098|MiY + O.23446IM2Y + 4.64976IM3Y, 
^Ats) = -0.20839(m?53)2 + 0.07828(m2^i)2 + 0.03346(m?^2)^ + 0.28668(771^3)^ + 0.28668(771^3)^ 



0^ 
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-0.17493(m^3)^ + 0.03346(m^3)2 + 0.00725|A°p - 0.01117|A°p - 0.18803|^°|^ 
+0.00192(A°A°* + v4°*v4°) - 0.00469(A°A°* + v4°*A°) + 0.07397(A°A°* + 
+0.00035(M°A°* + M°*A°) - 0.00932(M°A°* + M^* A^) + 0.01609(M°A°* + M^* A^) 
-0.00196(M°A°* + M°M°) - 0.01776(M2M^* + M°*Ag) + 0.03236(M°A°* + M°*A°) 
-0.02048(M3°A°* + M°M°) + 0.04148 (M3°A|^* + M^M^) - 0.04521(M3°A°* + M^M") 
-0.00595(M{'M°* + Mf M°) + 0.00716(M{'M3°* + Mf^M^) - 0.01505(M°M3°* + M^*Ml) 
+0.10639|MiY + 0.38267|M°P - I.72886IM3Y. (24) 

One can see from Eq. (24) that for At, Ah and S^, there are very large coefficients 
appearing before M3 or IM3 p which arise from the large gauge coupling of SU(3) group in the 
RGEs, Eqs. (4) and (5), and hence make them inevitably greatly depend on the gluino mass. 
If we consider the scenario with the mass spectrum below 1 TeV, we find that the effects of 
other high energy boundary values on the low energy spectrum are negligible as illustrated in 
Eq. (24) and actually it is gluino mass which mainly govern the SUSY mass spectrum (note 
that m° ~ or less than M° is required due to the constraint from cosmology [rH||) so that non- 
universality has a little influence on low energy mass spectrum. This is the screening effects 
induced by the large Yukawa couplings of top, bottom and tau for which a detailed discussion 



in the real parameter case has been given in Ref. |]13[. For At- and S^, the screening effects are 
weak since here M3 is not as important as for At, A^, and S;, and the coefficient of M2 can 
compete with that of M3 . In particular, for A^-, the coefficient of ^4° is the same order as that of 
M3 . That is, we are in the vicinity of the IRQFP in the example. Note that a crucial difference 
between our discussion and the one in Ref. [1^ is that we use the Yukawa couplings at mz scale 



induced from their corresponding current masses of quarks, while in Ref. [|T3l Yukawa couplings 



at the GUT scale much larger than the unified gauge couplings are assumed when the IRQFP 
behavior is discussed for large tan/?. Another difference is that we take the SUSY threshold 
at MsusY whereas in Ref. [|13| SUSY threshold effects was not considered. As a consequence 
we find that the value of tan/5 for which the IFQFP is reached depend weakly on the choice of 
MsusY- For example the theory exhibits IRQFP behavior at tan/5 ~ 59 for Msusy = 400 GeV 
and at tan/? ~ 60 for Msusv = 800 GeV. 

In Fig. 1(a) we illustrate how the dependences of At, A^ and Ar on their initial values, i.e., 
the coefficient clits) ( l=t,b,r) in Eq. (22), change with tan/?. One may find in the figure that 
the coefficient cKts) increases rapidly for tan/? less 5 and decrease slowly afterwards. This is 
because the top Yukawa coupling, yt, behaves as cx and for somewhat large tan/? it remains 
almost unchanged. So, according to the procedure described above, when tan/5 increases 
can not reach the value large enough to make At independent of ^4°. On the other hand, 
the coefficients cl and cj^. decrease considerably as tan/? increases in the whole range and drop 
sharply to zero when tan/5 approaches 59, which is exactly the case as expected since their 
corresponding Yukawa couplings behave as oc When tan/? reaches some large value ( ~ 59 
), Y^T- become large enough to make the dependences of Ab,r on their initial values, A^^, almost 
disappear, i.e., IRQFP is reached. In Fig. l(b-d) we plot the the dependences of T,t(ts), ^b(ts) 
and J^rits) on the initial values of their corresponding three parameters in the definitions of 
Tit, and S,-, i.e., fa{ts), k=t,b,r in Eq. (|23|) . One can find the similar behavior of the 
dependences on initial conditions when changing tan/5. It is clear from Fig. 1 that except for 
At and Sj, IRQFP behavior becomes more and more evident as tan/? approaches 59, i.e., the 
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point where Landau pole appear as we run bottom and tau Yukawa couplings up. 

In Fig. 2 we plot pi = Ai/M^^l = t,b and t) versus ag for tan/? = 58.7, running from the 
GUT scale down to Msusy scale as indicated by the axes of the strong coupling ag- One can 
clearly see the IRQFP behavior of piil = t,b and r) from the figure. In summary, the numerical 
results confirm that despite the CP violating phases introduced, the model still exhibit the 
IRQFP behavior for soft terms Ai and S^, as pointed out in section two. 

Now we come to a position to address the effects of CP violating phases on sparticle spec- 
trum. For the effects of phases of gaugino masses, and M2 , on the mass spectrum(because 
of R-symmetry, we take M3 real for convenience in the following discussion) , one can see from 
Eq. (24) that the only way for the phases of M° and M2 to be important is that their magni- 
tudes are much larger than M3 (say, at least one order of magnitude), which means a very heavy 
bino or chargino mass and so is not a preferable scenario. For the phases of trilinear terms A^, 
A'^ and A^, naively one may think that it is possible for them to play important role in the 
mass spectrum through the interference terms with M3 as illustrated in Eq. (24). However, 
notice that because the coefficients of the A'lM^* terms are much smaller than that of IM3P 
term ( e.g., for Et, the former are only 1/20 of the latter ) we have to choose |^°|(/ = t,b and r) 
much larger than M3 (say, at least one order of magnitude) in order to make the effects of their 
phases considerable so that one is led to far from IRQFP. Moreover, since trilinear soft SUSY 
breaking terms appear in the nondiagonal terms of squark or slepton mass matrices, the large 
Ai induced by very large ^41^ , together with terms proportional to /i tan/S, may result in so 
light stau that stau becomes LSP. Therefore, we can infer from the above discussion that the 
initial values of trilinear soft SUSY breaking terms cannot be much larger than the magnitudes 
of . Thus wc can conclude that in physically interesting regions of parameter space where 
the notion of IRQFP can be applied, CP violating phases of gaugino masses and soft trilinear 
SUSY breaking terms are not important for determining the mass spectrum through their roles 
played in RGEs. 

4 Conclusion 

In conclusion, we have shown in this paper another way to find the solutions of the RGEs of the 
soft SUSY breaking terms in the third generation sector of MSSM by directly solving the RGEs 
themselves. Compared with the method of the spurion expansion, our approach proves to be 
a simple and convenient way especially for models with complex soft trilinear SUSY breaking 
terms and gaugino masses. The results in the scenario with complex trilinear terms and gaugino 
masses turn out to support the notion of IRQFP. It follows that non-universality of trilinear 
couplings Ai and gaugino masses have no infiuence on sparticle spectrum at the IRQFP. We 
have studied the effects on mass spectrum of CP violating phases of trihnear terms and the 
non-universal gaugino masses and find that the effects on the mass spectrum are limited to be 
small in the physically interesting region of parameter space where the notion of IRQFP can 
be applied and physical requirements such as LSP being neutral particle are imposed. 
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Figure Caption 

Fig. 1(a) lines labelled by 1, 2 and 3 correspond respectively to coefficient of in At{ts), 
i.e., cl{ts) in (22), coefficient of in Ai,{ts), cl{ts), and coefficient of A'^ in Ar{ts), cl{ts), as 
functions of tan/3. In 1(b) lines 1 and 2 refer respectively to coefficients of (m^g)^ and (m^g)^ 
in Y:t{ts)- In 1(c) lines labeled by 1, 2 and 3 are respectively coefficients of (rngg)^, {mPjj-^)'^ and 
(m^g)^ in S6(t5). In 1(d) lines 1 and 2 correspond respectively to coefficients of (m^^^)^ and 
{m%f in E.(t5). 

Fig. 2 Three dimensional plots of pi = Ai/M^{l = t,b and r) as functions of as for tan/3 = 
58.7. The other parameters are chosen as M° = M° = 500 GeV, = 400 GeV and Mq = 800 
GeV. Initial values of |p/|'s are chosen as 1 or 2 and the phases of p°'s are taken asiX7r/4(i = 
0,...,7). 
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